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On different cascade-speeds for longitudinal and transverse velocity increments 
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We address the problem of differences between longitudinal and transverse velocity increments in 
isotropic small scale turbulence. The relationship of these two quantities is analyzed experimentally 
by means of stochastic Markovian processes leading to a phenomenological Fokker- Planck equa- 
tion from which a generalization of the Karman equation is derived. From these results, a simple 
relationship between longitudinal and transverse structure functions is found which explains the 
difference in the scaling properties of these two structure functions. 
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Substantial details of the complex statistical behaviour 
of fully developed turbulent flows are still unknown, 
cf- 1^ Q|' One important task is to under- 

stand intermittency, i.e. finding unexpected frequent 
occurences of large fluctuations of the local velocity 
on small length scales. In the last years, the differ- 
ences of velocity fluctuations in different spatial direc- 
tions have attracted considerable attention as a main 
issue of the problem of small scale turbulence, see for 
example HHSIllillSElllllIl- For local 
isotropic turbulence, the statistics of velocity increments 
[v(x + r) — v(x)] e as a function of the length scale r is 
of interest. Here, e denotes a unit vector. We denote 
with u{r) the longitudinal increments (e is parallel to r) 
and with v{r) transverse increments (e is orthogonal to 

In a first step, this statistics is commonly investigated 
by means of its moments {u'^{r)) or (v"(r)), the so-called 
velocity structure functions. Different theories and mod- 
els try to explain the shape of the structure functions 
cf. [3- Most of the works examine the scaling of the 
structure function, (u") oc r^' , and try to explain in- 
termittency, expressed by — n/3 the deviation from 
Kolmogorov theory of 1941 For the correspond- 

ing transverse quantity we write (v") oc r^' . There is 
strong evidence that there are fundamental differences 
in the statistics of the longitudinal increments u{r) and 
transverse increments v{r). Whereas there were some 
contradictions initially, there is evidence now that the 
transverse scaling shows stronger intermittency even for 
high Reynolds numbers 0, 0| . 

A basic equation which relates both quantities is de- 
rived by Karman and Howarth |l7j . Assuming incom- 
pressibilty and isotropy, the so called first Karman equa- 
tion is obtained: 
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(1) 



Relations between structure functions become more and 



more complicated with higher order, including also pres- 
sure terms 0, 0| . 

In this paper, we focus on a different approach to char- 
acterize spatial multipoint correlations via multi-scale 
statistics. Recently it has been shown that it is pos- 
sible to get access to the joint probability distribution 
p{u{ri),u(r2)-, ■ ■ ■ ,u(rn)) via a Fokker-Planck equation, 
which can be estimated directly from measured data 
[lllHiJ. For a detailed presentation see [22]. This 
method is definitely more general than the above men- 
tioned analysis by structure functions, which characterize 
only the simple scale statistics p{u{r)) or p{v{r)). The 
Fokker-Planck method has attracted interest and was 
applied to different problems of the complexity of tur- 
bulence like energy dissipation [231 |2al, universality 
turbulence ^ and others [H lU llj ill ^3, EH ■ The 
Fokker-Planck equation (here written for vector quanti- 
ties) reads as 
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A, p(u,r|uo,ro). 
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{i denotes the spatial component of u, we fix i = 1 
for the longitudinal and i = 2 for the transverse incre- 
ments.) This representation of a stochastic process is 
different from the usual one: instead of the time t, the 
independent variable is the scale variable r. The minus 
sign appears from the development of the probability dis- 
tribution from large to small scales. In this sense, this 
Fokker- Planck equation may be considered as an equa- 
tion for the dynamics of the cascade, which describes 
how the increments evolve from large to small scales un- 
der the influence of deterministic {D^^") ) and noisy (Z?^^)) 
forces. The whole equation is multiplied without loss 
of generality by r to get power laws for the moments 
in a more simple way, see also Eq. (Q. Both coeffi- 
cients, the so-called drift term D^^^\u,r) and diffusion 
term Dlj{u,r), can be estimated directly from the mea- 
sured data using its mathematical definition, see Kol- 
mogorov 1931 [pl and [2l|13,Ell2l- With the notation 
AUi{r, Ar) := Ui{r — Ar) — Ui{r) the definitions read as: 
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Dl'\u,r) = ^lim^^(Ai7,(r, Ar))|uM=u(r), 
Difin,r)= (3) 

Here we extend the analysis to a two dimensional 
Markov process, relating the longitudinal and trans- 
verse velocity increments to each other. The resulting 
Fokker-Planck equation describes the joint probability 
distribution p{u{ri), v{ri); . . . ; M(r„), w(r„)). Knowing 
the Fokker-Planck equation, hierarchical equations for 
any structure function {f{u, v)) = J f{u, v)p{u, v, r)dudv 
can be derived by integrating Eq. 

^ (4) 

+ m{u'"-'^v"D[^^ (u, V, r)) + n(u"'v''-'^ D^^^ {u, v, r)) 

+ mn{u"'-^v"~^D^^^\u,v,r)), 

which we take as a generalization of the Karman equa- 
tion. 

Next, we describe the experiment used for the subse- 
quent analysis. The data set consists of 1.25-10^ samples 
of the local velocity measured in the wake behind a cylin- 
der (cross section: D=20mm) at a Reynolds' number of 
13236 and a Taylor- based Reynolds' number of 180. The 
measurement was done with a X-hotwire placed 60 cylin- 
der diameters behind the cylinder. The component u 
is measured along the mean flow direction, the compo- 
nent V transverse to the mean flow is orthogonal to the 
cylinder axis. We use Taylor's hypothesis of frozen tur- 
bulence to convert time lags into spatial displacements. 
With the sampling frequency of 25kHz and a mean ve- 
locity of 9.1 m/s, the spatial resolution doesn't resolve 
the Kolmogorov length but the Taylor length A = 4.85 
mm. The integral length is L =137 mm. 

From these experimental data the drift and diffusion 
coefficients are estimated according to eqs. Q as de- 
scribed in [33, see also [s^. As an example, the diffusion 
coefficient D^^ {u,v,r = L/A) is shown in Fig. ^ To use 
the results in an analytical way, the drift and diffusion 
coefficient can be well approximated by the following low 
dimensional polynoms, which will be verified by recon- 
structed structure functions as it is shown below: 

D[^\u,v,r) = d'l{r)u (5) 

Di'\u,v,r) = dUr)v 

D^^^iu,v,r) = dn{r)+d'l,{r)u + d^nrW +dri{ry 

D'i^iu, V, r) = d22{r) + d!^^{r)u + dl^^{r)u^ + dll{r)v^ 

D^^^{u,v,r) = di2{r)+d'^2{r)u + d'i:^{r)uv. 




FIG. 1: Diffusion coefficient (u,v,r) for r = L/4 esti- 
mated from experimental data. Note that the quadratic con- 
tributions are responsible for intermittency. The asymmetry 
in It-direction is related to the non- vanishing skewness of the 
probability distribution of longitudinal velocity increments. 



In order to show that the Fokker-Planck equation with 
these drift and diffusion coefficients can well characterize 
the increment's statistics, one has to verify that the evo- 
lution process of u{r) and v(r) is a Markov process and 
that white and Gaussian distributed noise is involved. 
The Markov property can be tested directly via its def- 
inition by using conditional probability densities [23| or 
by looking at the correlation of the noise of the Langevin 
equation |22 . For our case we have verified that the one- 
dimensional processes of the longitudinal and transverse 
increments are Markovian |4ll | , thus the two- dimensional 
processes should be Markovian, too |3^. As an alterna- 
tive approach to verify the validity of the Fokker-Planck 
equation we have solved numerically the hierarchical eq. 
Q for {u^) using the above mentioned coefficients. To 
close the equation, we have used the moment (w^^^u^) 
from the experimental data. In Fig. |21 the integrated 
longitudinal structure functions are given in comparison 
with the structure functions directly calculated from the 
data (with n =2, 4, 6 and m =0). These results we take 
as the evidence that the Fokker-Planck equation charac- 
terizes the data well and can be used for further inter- 
pretations. 

For the drift coefficient, which is the deterministic part 
of the cascade dynamics, the process decouples for the 
different directions. The drift and diffusion coefficients 
are symmetric with respect io v —v. Remark, in con- 
trast to the statistics of the longitudinal increments, the 
transverse one is symmetric and show for example no 
skewness (u'^) = 0. Furthermore, quadratic terms oc- 
cur in the diffusion coefficients. Intermittency results 
from the quadratic terms cZ"" and (^221 a-U other terms 
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FIG. 2: Even longitudinal structure functions up to order 6 
calculated directly from data (symbols) are compared to the 
results obtained from numerical solutions of Eq. using 
the experimentally estimated coefficients of the Fokker-Planck 
equation (solid lines; see also Eq. @ and Fig. . 





1 1 1 

i/22(2/3r) 


1 1 






















f/7(2/3r) 
-i , 1 


"*" + + 

1 



0.4 



0.3 



0.0 



(b) 



rf„""(r) 
i/22"'(2/3r) 



01 - I^**.w.V*W.>**"'*^-* 



0.00 0.05 0.10 0.15 0.00 0.05 0.10 0.15 

r [m], Wr r [m], 2/3r 

FIG. 3: The expansion coefficients of the drift and diffusion 
coefficients in dependence of the scale r. The abscissa are 
rescaled for the transverse coefficients. The corresponding 
longitudinal and transverse coefficients coincide which each 
other a) apart from the intermittency terms b). 
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FIG. 4: Longitudinal (symbols) and corresponding transverse 
(lines) structure functions. The arrows point from the struc- 
ture function {«"(?■)) to the corresponding structure function 
(u"(2/3 r)) with contracted abscissa. The transverse struc- 
ture functions with a contracted abscissa are close to the lon- 
gitudinal ones. 



Next, let us suppose that the structure functions scale 



with a power law, 



(r)) 



c"r^" and 



even though our measured structure functions are still 
far away from showing an ideal scaling behavior |2(i| |. 
With exemption of the differences between di"(r) and 
^22(^)1 can relate the structure functions according 
to the above mentioned rescaling: (w"(r)) = (w"(|r)) = 



t ' ■ — ^rd'') ' • We end up with the relation = ^" 



and % = (I 

c, V 2 



Note that the c" constants are related 



to the Kolmogorov constants. For n = 2 and 4 we obtain 
Cf/cf « 1.33 and cj/cf « 1.72, which deviates less than 
3% from the value of c^/cf = 4/3 and c^/cf = 16/9 given 



act against intermittency. 

The r-dependence of related longitudinal and trans- 
verse d-coefficients {d" and d!^ etc.) coincides if the 
abscissa are rescaled: diongC'') ~ c?transv(|''), see Fig. 
13 The only exception is the coefficient (^22(^)1 whereas 
d1^{r) « d^l'(r) « ^22 (^) ~ const.. We interpret this 
phenomenon as a faster cascade for the transverse in- 
crements. It can be seen from the hierarchical equation 
that this property goes over into structure functions 
of arbitrary even order, (t;"(r)) w (M"(|r)). Only the 
small coefficients dii and ^22 break this symmetry, be- 
cause they belong to different odd, and therefore small, 
moments. In Fig. ^ the structure functions of order 
2, 4 and 6 are plotted with respect to this rescaled 
length The structure functions are normalized by 
(a")* = (q!")-^t^^|^, with a either u oi v. 

The observed phenomena are consistent with the 
Karman equation J^l, if the Karman equation is inter- 
preted as a Taylor expansion 



(6) 



{uHr + lr))^{uH^r)). (7) 



At last we discuss the use of ESS (extended self- 
similarity |37|l Vwith respect to transverse velocity com- 
ponents. In [81 13^13^ the authors plot (w") against {\u^\) 
and obtain that the transverse exponents is smaller than 
the longitudinal one, < Ci"- In Fig- Elthe fourth struc- 
ture functions are plotted against which shows 
clearly that < If the transverse structure function 
is plotted as a function of (|u'^(|r)|), this discrepancy 
vanishes. Notice that these properties are due to a none 
existing scaling behavior. It is evident that our rescal- 
ing does not change the exponents in case of pure scaling 
behavior. 

To conclude the paper, we have presented experimental 
evidence that the statistics of longitudinal and transverse 
increments is dominated by a difference in the "speed of 
cascade" expressed by its r dependence. Rescaling the 
r dependence of the transverse increments by a factor | 
fades the main differences away. Thus the longitudinal 
and transverse structure functions up to order 6 coincide 
well. A closer look at the coefficients of the stochastic 
process estimated from our data shows that the multi- 
plicative noise term for the transverse increments ^22 ^'^d 
the symmetry breaking terms ^-nd ^22 do not follow 
this rescaling. These coefficients may be a source of dif- 
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ferences for the two directions, but for our data analyzed 
here this effect is still very small. 



Finally, we could show that our findings on the rescal- 
ing are consistent with the Karman equation and that 
longitudinal and transverse Kolmogorov constants of the 
structure functions up to order four can be related con- 
sistently with our results. 



FIG. 5: The fourth order longitudinal (x) and transverse (+) 
structure function is plotted against (liip), the third moment 
of the longitudinal increments' modulus (ESS plot). If the 
abscissa of the transverse structure function is rescaled by a 
factor of 2/3, both curves fall one upon the other (squares). 
The differences between the exponents also vanish. 



We acknowledge fruitful discussions with R. Friedrich 
and teamwork with M. Karth. This work was supported 
by the DFG-grant Pe 478/9. 



[1] A. S. Monin and A. M. Yaglom, Statistical Fluid Mechan- 
ics: Mechanics of Turbulence, vol. 2 (The MIT Press, 
1975). 

[2] U. Frisch, Turbulence (Cambridge University Press, 
1995). 

[3] A. Tsinober, An Informal Introduction to Tubrulence 

(Kluwer Academic Publ., 2001). 
[4] M. Nelkin, Am. J. Phys. 68, 310 (2000). 
[5] R. Antonia, B. Pearson, and T. Zhou, Phys. Fluids 12, 

3000 (2000). 

[6] I. Arad, B. Dhruva, S. Kurien, V. L'vov, I. Procaccia, 

and K. Sreenivasan, Phys. Rev. Lett. 81, 5330 (1998). 
[7] S. Chen, K. Sreenivasan, M. Nelkin, and N. Cao, Phys. 

Rev. Lett. 79, 2253 (1997). 
[8] B. Dhruva, Y. Tsuji, and K. Sreenivasan, Phys. Rev. E 

56, R4928 (1997). 
[9] T. Gotoh, D. Fukayama, and T. Nakano, Phys. Fluids 

14, 1065 (2002). 
[10] S. Grossmann, D. Lohse, and A. Reeh, Phys. Fluids 9, 

3817 (1997). 

[11] J. Laval, B. Dubrulle, and S. Nazarenko, Phys. Fluids 

13, 1995 (2001). 
[12] W. van de Water and J. Herweijer, J. Fluid Mech. 387, 

3 (1999). 

[13] R. Hill, J. Fluid Mech. 434, 379 (2001). 
[14] A. N. Kolmogorov, Dokl. Akad. Nauk SSSR 32, 16 
(1941). 

[15] A. Kolmogorov, J. Fluid Mech. 13, 82 (1962). 
[16] X. Shen and Z. Warhaft, Phys. Fluids 14, 370 (2002). 
[17] T. von Karman and L. Howarth, Proc. Roy. Soc. A 164, 
192 (1938). 

[18] R. Hill and O. Boratav, Phys. Fluids 13, 276 (2001). 
[19] R. Friedrich and J. Peinke, Physica D 102, 147 (1997). 
[20] R. Friedrich and J. Peinke, Phys. Rev. Lett. 78, 863 
(1997). 

[21] P. Marcq and A. Naert, Phys. Fluids 13, 2590 (2001). 
[22] C. Renner, J. Peinke, and R. Friedrich, J. Fluid Mech. 

433, 383 (2001). 
[23] A. Naert, R. Friedrich, and J. Peinke, Phys. Rev. E 56, 



6719 (1997). 

[24] C. Renner, J. Peinke, R. Friedrich, O. Chanal, and 

B. Chabaud, arXiv: physics/0211121. 
[25] P. Marcq and A. Naert, Physica D 124, 368 (1998). 
[26] C. Renner, J. Peinke, R. Friedrich, O. Chanal, and 

B. Chabaud, Phys. Rev. Lett. 89, art. no. 124502 (2002). 
[27] J. Davoudi and M. Tabar, Phys. Rev. Lett. 82, 1680 

(1999). 

[28] T. Frank, Physica A 320, 204 (2003). 

[29] I. Hosokawa, Phys. Rev. E 65, art. no. 027301 (2002). 

[30] M. Ragwitz and H. Kantz, Phys. Rev. Lett. 8725, art. 

no. 254501 (2001). 
[31] F. Schmitt and D. Marsan, Eur. Phys. J. B 20, 3 (2001). 
[32] A. N. Kolmogorov, Math. Ann. 104, 415 (1931). 
[33] H. Risken, The Fokker-Planck Equation (Springer- 

Verlag, 1989). 

[34] R. Friedrich, C. Renner, M. Siefert, and J. Peinke, Phys. 
Rev. Lett. 89, art. no. 149401 (2002). 

[35] J. Peinke, C. Renner, M. Karth, M. Siefert, and 
R. Friedrich, in Advances m Turbulence IX, Proceedings 
of the Ninth European Turbulence Conference, edited by 
I. Castro, P. Hancock, and T. G. Thomas (CIMNE, 
2002), p. 319. 

[36] R. Antonia, M. Ould-Rouis, Y. Zhu, and F. Anselmet, 

Europhys. Lett. 37, 85 (1997). 
[37] R. Benzi, S. Ciliberto, R. Tripiccione, C. Baudet, F. Mas- 

saioh, and S. Sued, Phys. Rev. E 48, R29 (1993). 
[38] B. Pearson and R. Antonia, J. Fluid Mech. 444, 343 

(2001). 

[39] R. Antonia and B. Pearson, Europhys. Lett. 40, 123 
(1997). 

[40] Note we focus on transverse increments for which r is in 
direction to the mean flow and v is perpendicular to the 
mean flow. An other possibility is to chose r perpendicu- 
lar to the mean flow, which possibly results in a different 
behavior in case of anisotropy |12.ll6| . 

[41] to be published 



